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Abstract 

We solve the Klein-Gordon and Dirac equations in an open cosmological 
universe with a partially horn topology in the presence of a time dependent 
magnetic field. Since the exact solution cannot be obtained explicitly for 
arbitrary time-dependence of the field, we discuss the asymptotic behavior of 
the solutions with the help of the relativistic Hamilton-Jacobi equation. 
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I. INTRODUCTION 



During the last years a large amount of observational data has been reported showing 
that our universe is almost isotropic and homogeneous. The study of the structure of the 
Cosmic Microwave Radiation leads us to conclude that the ratio of the total density to the 
critical density of the universe flo is likely to be close to one 0-^], favoring a spatially flat 
Robert son- Walker metric over other topologies. 

It is well known that general relativity is a local metrical theory and therefore the cor- 
responding Einstein field equations do not fix the global topology of spacetime and con- 
sequently the universe may have compact spatial sections with a nontrivial topology ||||, 
then the observational data does not rule out the possibility that our universe possesses a 
hyperbolic topology 

The study of cosmological models with nonstandard topologies is not new and goes back 
to the works by Zelmanov PJTOl, showing that upon different coordinate transformations, 



spatially closed or flat sections can be transformed into hyperbolic sections and vice versa. 
The line element associated with an spatially open Friedman universe has the form 

ds 2 = a 2 (r]) [-drf + dr 2 + sinh 2 (d6 2 + sin 2 6d(f) 2 )] (1) 

Making the coordinate transformation 

e~ z = cosh r — sinh r cos 6, e~ z x = sin 9 cos sinh r, e~ z y = sin 6 sin <p sinh r (2) 

the metric ([]]) becomes 

a~ 2 (ri)ds 2 = -d V 2 + dz 2 + e~ 2z (dx 2 + dy 2 ) (3) 

The topology is induced by identifying points periodically along x and y by (x, y) = (x + 
b,y + h), where b and h are constant, to create a two dimensional torus. The torus is 
stretched by the factor e~ 2z along the z axis to create a toroidal horn. The comoving proper 
area of the torus is e~ 2z bh and depends on location along the z axis. The global topology 
induces global inhomogeneity as well as global anisotropy [H . 



The study of quantum effects in cosmological universes with a nontrivial topology allows 
us a deeper understanding of the properties of different scenarios and which of them can 
describe our universe. In this direction the metric (|3]) represents a very interesting scenario 
in order to discuss particle production and propagation of perturbations in cosmology. 

After the publication of the pioneering article by Schrodinger discussing particle 
production in a deSitter universe, many articles have been published on the problem of 
quantum effects in cosmological scenarios [|T2|— ^L4 1 , most of them dealing with a Robertson- 
Walker line element with spatially flat topology. This particularly simple line element, which 
is the most used in inflationary models, permits one to compute the Green function as well 



as to solve the relativistic wave equations [15 17 



In order to study quantum processes in curved space-times one has to fulfill a preliminary 
step which consists in having a description of the single-mode solution of the relativistic 
particles or perturbations in those background fields, i.e., exact solution of the relativistic 
scalar and spinor wave equations. In the literature we have at our disposal different methods 
for solving relativistic wave equations in curved spaces and in curvilinear coordinates; among 
them the method of separation of variables is one of the most widely used. [ITSH2TI . 



It is the purpose of the present article to solve the Klein-Gordon and Dirac equations in 
the Friedman universe associated with the metric (f|) in the presence of a time dependent 
magnetic field. In order to solve the Dirac equation we apply the algebraic method of 



separation of variables p0|-p4|. We compare the solutions with those of obtained after solving 
the relativistic Hamilton Jacobi equation. The article is structured as follows: In Sec. II 
we solve the relativistic Hamilton- Jacobi equation in an open cosmological universe with a 
horn topology. In Sec. Ill we separate variables and solve the Klein-Gordon equation. In 
Sec IV, using the algebraic method of separation of variables, we reduce the Dirac equation 
to a system of first order coupled differential equations that we solve in terms of special 
functions. Finally, in Sec. V we briefly discuss the results reported in this article. 
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II. SOLUTION OF THE HAMILTON- J ACOBI EQUATION 



The covariant generalization of the Hamilton- Jacobi equation has the form [25 

dS 



9 



eA 



M z = 0, 



(4) 



dx a "J ydx 13 

where the contravariant metric, A a is the vector potential and M is the mass of the 

particle. Here and elsewhere we adopt the convention c = H = 1. 

Let us introduce an electromagnetic field associated with the vector potential 



(5) 



where the index fi = is associated with the evolution parameter r\ and fi = 1,2,3 correspond 
to the space coordinates x, y, z respectively. Looking at the relativistic invariants 



Az 



2 ^ 



B — E 



dy 



f^f; v = o, 



(6) 



(7) 



and taking into account that only is different from zero, we notice that ([5]) corresponds 
to a non constant magnetic field B, directed along the z axis, with strength. 

1 dA 1 (y) 



B 



dy 



(8) 



whose value is inversely proportional to the expansion factor a(r)) 2 . 

The line element (||) is a Stackel space [26|, and the Hamilton- Jacobi equation (|J) is 
completely separable in (|3|) in the presence of the vector potential (|5]), therefore we can look 
for a solution in the form 



S = k x + S y (y) + S z {z) + S v {r)). 



(9) 



Substituting ([5]) into Eq. 

(h-A^y)) 2 , 



we obtain 



-2z 



-2z 



dy 



+ 



dSz 
dz 



dS'q 

drj 



- M 2 a(r]) 2 = 0. 



(10) 
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Equation fllCf ) reduces to the following system of differential equations: 



§v + ^=4 (id 



dS x 2 



di] 



(12) 



^-My)) 2 +[^-) =k 2 xy , (is) 

where fc^, and k\ are separation constants. 

In the absence of electromagnetic interaction, we have that Ai(y) = and the solution 
of Eq. (fl3|) takes the form 



S y = ± 1 Jk* u -k*y = ±k y y. (14) 

where we have introduced the constant k y . Equation (|14]) can also be derived looking at the 
symmetry between the torus coordinates x and y in the line element (||]) and Eq. (|lOj) when 
A 1 (y) = 0. 

When the vector potential has the simple form A\(y) = Aiy, the magnetic field reads 

2z 

B = Q e ^ 2 \Ai\ and the function S y (y) is 

S y (y) = -^^Jk* y -(k x -A iy y + ^ a rct a n A vV~ k * =i (15) 

The solution of Eq. ( |TTD can be expressed in terms of elementary functions as follows: 



-i k2 z-kl y exp(2z) 



S z = yjkl - kl y exp(2^) - k z tanh" 1 \ / ^ ^ " v ' . ( I (i ) 

The solution of Eq. ( |i~2"D can be written as 

S V = ±J Vk 2 z ~ M 2 a{ri) 2 dr], (17) 

whose explicit form in terms of elementary functions will depend on a particular choice of 
the expansion function ot{rj). 
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Since we have been able to solve the Hamilton- Jacobi equation in the Stackel space given 
by we can construct the quasiclassical modes of the relativistic wave equations through 
the identification 



_> e iS — e ±iJ\/k%-M 2 a(r)) 2 d,ri e ik x x e iSy e iS z ^ ^g-j 

where S z and S y take the following values at the asymptotes 

*Sz(oo) > ik xy e , S^—oo) >■ fc^-z, (19) 

<Voo) T* ' ^ 

When the electromagnetic interaction is not present we have that S y = exp(ik y y). 

III. SOLUTION OF THE KLEIN GORDON EQUATION 

The covariant generalization of the Klein Gordon equation in curved space-time has the 
form [T2\ 



g afi (V Q - ieA a ) (Vp - ieA p ) $ — (M 2 + £R)$ = 0, (21) 

where 

V Q is the covariant derivative, R is the curvature scalar and £ is a scalar dimensionless 
coupling constant which takes the value £ = 1/6 in the conformal case and £ = when a 
minimal coupling is considered. The value of the R for the metric (^) is 

—a(<n) + d2a(??) 

* = 6 air,)^ - < 22 > 
Substituting the metric associated with the line element (|j) into the Klein Gordon equation 



(21]) one obtains 



2 ,<9 2 $ 22 <9 2 $ 0$ r3 2 $ 3$ da 1 224 , n . 2z d§ . . . „ n 

e 7TT + e 7TT - 2 IT + 7TT ~ 2 7T1 3 ~ e Mv) ^ ~ ^e 2z —A x {y) - M 2 a 2 $ = 0, 

ox z Oy z oz Oz z or] arj cr ox 

(23) 



where we have chosen to work with a minimal coupling £ = 0. The Klein Gordon equation 



is completely separable in (EI), therefore we look for its solution in the form. 



$ = H(rj)Z(z)Y(y)e 



(24) 



Substituting ( |2"4"D into Eq. (^Tj) we reduce the problem of solving the Klein-Gordon equation 
to that of finding solutions of the following set of ordinary differential equations 

d 2 Y 



dy 2 



-((k x -A 1 (y)) 2 -k*)Y = 0, 



(25) 



d 2 Z dZ 



dz 2 dz 



2— - (A + k 2 e 2z )Z = 



(26) 



d 2 H dH da , , , „ , Q , „ 

with A 2 and k 2 as separation constants. For A(y) = A\y the solution of Eq. 
expressed in terms of Whittaker functions [571 as follows 



(27) 



']) can be 



(28) 



where 



Aiy - k z 



(29) 



and C\ and are arbitrary constants. In the absence of electromagnetic field the solution 
of Eq. fl25|) reduces to 



(30) 



The solution of Eq. © is J28 



Z = C 3 e*H^(ike*) + C 4 e* H%^(ike*) , (31) 



r(2) 



where Hj) 1 \z) and Hl 2 \z) are the Hankel functions and C3 and C4 are arbitrary constants. 
We can also express the solution of (|26|) in terms of Bessel functions J v {z) as 



Z = D z e z J VT ^{ike z ) + D 4 e z J_ 



[ike z 



(32) 



where D 3 and D± are arbitrary constants. 
After introducing the function ^.(77): 



a{rj)' 



(33) 



Eq. 



reduces to 



d 2 h 



+ (a^M 2 - A 2 - t^L)h = 0. 



(34) 



dr/ 2 drf 

In order to analyze the asymptotic behavior of the solutions of the Klein-Gordon equation 
(PH) we make use of the asymptotic behavior of the Hankel functions [27 



1/2 



= j(2-7T^/2-7r/4) 



1/2 



-i(z— -nv/2— 7r/4) 



(35) 



as 2 — ► 00, and the behavior of J v (z) as 2 — > |^| 



r i + 



(36) 



The asymptotic behavior of the Whittaker function Wk,p(z) for large values of z is [^8 



W^z) - e- z ' 2 z\ 



(37) 



and the function M^iz) has the following asymptotic behavior as z — > 



AM* 



(38) 



An approximate solution of Eq. ( p4|) can be obtained provided that the expansion parameter 
ol{j]) satisfies the conditions of validity of the adiabatic approximation. In this case one has 
that h{rj) has the form 



1 Pi, 

==exp(±i / W{rj)drj 



(39) 



with 



W{r 1 ) 2 = u{r 1 ) 2 [l + 5 2 (r))u- 2 + ...] 
8 



(40) 



where the function u {rj) has the form 

co( V y = a 2 (r))M 2 - A 2 - (41) 
5 n (r]) is a function of uj(r]) and its derivatives at i] up through u^(rj) and # n (77) is bounded 



as lu — > oo. The solution of the Klein-Gordon equation (EI) can be written as 



exp(±z f Ja 2 (ri)M 2 - A 2 - ^#^77 
$ = J V ^ ^-Z(*)Y(y)e^. (42) 



Let us analyze the asymptotic behavior of (fi^) as 7/ — > 00 and z — > —00. Using (|i~9|) and ( |36"D 
we obtain that, when the electromagnetic interaction is switched off, the mode solutions of 
Eq. (pip take the asymptotic form 



eX P (±7 r A / a 2( r/ ) M 2_ A 2_rf^?) d 

V2a{ V ) 3 / 2 ' 

Analogously, we have that in the presence of the electromagnetic potential the mode solutions 
of Eq. (|2l|) take the following asymptotic form 



exp(±z ( Ja 2 (rj)M 2 - A 2 - ^^dr] , 

For large positive values of z we have that the asymptotic behavior of $ is 



exp(±z ( Ja 2 (r])M 2 - A 2 - ^^cfy 

v/2a(r/) 3 / 2 V ; 

From Eq. (H4T) we can identify the quasiclassical modes as y — > oo and 2 ^ -00 as 



^d«.(*— 00) = ^fJ±VT+^^ keZ ) v ~ ll2w J^± (v 2 )e ik * x . (46) 



Analogously, from Eq. (f45|) we have that the quasiclassical modes as y — ► 00 and z — > 00 
are 

$da SS (^oo) = ^e z H { ^^(ike z )v-V 2 Wj^Jv 2 )e ik * x . (47) 
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IV. SOLUTION OF THE DIRAC EQUATION 



The Dirac equation is a system of coupled partial differential equations which is separable 
in a very restricted set of metrics. Among the spacetimes where the separability of the 
Klein-Gordon and Dirac equations has been studied one can mention the Stackel spaces 
|26|j , which are those metrics where the Hamilton- Jacobi equation is separable. Nevertheless 



recently it has been shown that this condition is neither necessary nor sufficient in order 



to guarantee a complete separability of variables in the Dirac equation (see Ref. [3^] and 



references therein). A systematic classification of the gravitational backgrounds where the 



Dirac equation is separable with the help of the algebraic method is presented in ref. pO 
The line element (|3]) belongs to this family and consequently one can apply the algebraic 
method of separation. 



The covariant generalization of the Dirac equation in curved space-time is |12| , |33 

T (d a -T a - ieA a ) § + M§ = 0, (48) 
where the curved Dirac matrices 7" satisfy the commutation relation 

{r^}=2g^, (49) 
and r a are the spin connections p3l 



r 1 

1 a — Tfl'/LtA 



dx" P ™ 



s^, (50) 



where 



s ^ = I { yt - rr) , (51) 

and the matrices 6", establish the connection between the Dirac matrices 7^ on a curved 
space-time and the flat Dirac matrices 7^ as follows: 

% = b a ^ a , r = a^- (52) 

Since the line element (^) is associated with a diagonal metric, we can work in the diagonal 
tetrad gauge for 7^: 
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7° 



7 - , f 



7 a 



r 



-,7 3 



7 



0(77)' ' a(rj)e~ z ' ' a(?7)e~ z ' ' 0(77) ' 
substituting (Q) into ( |50"D we obtain that the spinor connections are 

1 e~ ; 



2a( V ) 



(53) 



(54) 



1 e 



2a( V ) 



(55) 



1 docjrj) 1 34 

is = : ^T7 7 ) 1 4 



0. 



2 G?77 «(??) 

Substituting (K3)-(E6]) into (H3) we find that the Dirac equation takes the simple form 



(56) 



K + 7V (I " • 4l(! ' ) ) + ^fy + ^1 + Ma{ "> > * = 



where we have introduced the spinor \l/ 



(57) 



# = a(r/)~ 3/ V^. 



(58) 



Regarding Eq. ( 57Q we should mention that it does exhibits a nonfactorizable structure 
22| , |3"4]| . In order to solve Eq. fl5?D we apply the algebraic method of separation of variables 
. The method consists in rewriting the Dirac equation ( |5TD as a sum of two first order 



(59) 



differential operators Ki,K2 satisfying the relation 



with 



{). {/m • /v 2 }<I 



(60) 



and 



K x {x, y)$ = <i 7 2 ^ + 7 1 (Sj; - iMv) ) \ 7 3 7°$ = ik$, 



(61) 



#2(3, 77)$ = e z { 7°^ + 7 3 ^ + Ma(r)) \ 7 3 7°$ = -ik$. 



(62) 



11 



It should be noticed that using the pairwise scheme of separation one has been able to reduce 
the problem of solving the Dirac equation to finding solutions of the decoupled system of 
Eqs. (^Jl]) and (0). A further problem arises when we try to separate variables in Eq. (|62|). 
Here it is not possible to reduce the problem to a set of two commuting first order differential 
operators. In order to separate variables in Eq. flB2|) we re-write it in the following form: 

mm 



L l7 V + L 2 )$ = 0, (63) 



where Li and L 2 are two commuting differential operators given by the expressions 

Li — 7° + Ma(r]), (64) 

L 2 = ^^-+ike\ (65) 



In order to separate variables in Eq. fl63|) we introduce the auxiliary spinor 3^ 

L l7 3 7 ° + L 2 )y = $, (66) 



where the differential operator L 2 is given by the expression 

L 2 = 1 °-^--tke z . (67) 
Substituting (EB) into ( |5"3"D we obtain that y satisfies the following equation 

Mi + M 2 }y = 0, (68) 



with 



M ly M : 



0, and 

d 2 



Mx + \)y= (^-^ - i^ke z + k 2 e 2z + Xj y = 0, (69) 

(M 2 -X)y=(^ + 7 °M^1 + M 2 a 2 ( V ) - A) y = 0, (70) 

where A is a separation constant. Introducing the new variable u = 2ke z , we have that Eq. 
(|69|) can be written as 

12 



d 2 i o 1 .1 

1 7 h (- 

du 2 2u 1 4 l 4 



A)- 



5 = 0, 



where 



u 



- x ' 2 s = y. 



(71) 



(72) 



Choosing the following representation of the Dirac matrices |36 

\ 



7 



V° V 



ct j 

(T- 7 ' 



1 < J < 3 



(73) 



we readily obtain that the spinor $ has the following structure 

d 



ctit^ itr^kx - A x {y)) 

dy 



d 



$ 2 = ik&i, 



0(2/) 

i^ 3 <Ky) 



\ 



exp(ik x x) 



(74) 



(75) 



(76) 



/ 



where 



(77) 



y B(y) j 

Using the representation (^) we obtain that the solution of Eq. ([711) can be written in 
terms of Whittaker functions 

Si,2 = D 1 W_ 1/2> ^(u) + D 2 M_ 1/2 ~(u), S 3A = D 3 W 1/2 ^(u) + D 4 M 1/2 ^(u) (78) 



where D\, D 2 , D 3 , D± do not depend on the variable u. Looking at ( 170 ) and ( J7T] ) we have 
that, for regular solutions at u — 0, the spinor 3^ has the following structure: 

( a{y)cMu- l/2 M +hV - x {u)\ 
b(y)cMu~ l,2 M +lv - x (u) 
c{y)c 2 {ri)u- l l 2 M_ lvl {u) 

\ d(y)c 2 (7 7 ) W - 1 /2M_iy X (w) y 



3^ 



exp(z/c :r x) 



(79) 
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Substituting (|79|) into (|66D and noticing that Eq. ([7(]) is equivalent to the following system 
of equations 



d_ 

dr] 



- iMa(r))) d(r)) = Vlc 2 (v), 



(80) 



-|- + iMa(ri)^J c 2 (r}) = ^/Ic^r)), 
we obtain that the spinor $ has the following structure 



^ A(v) Cl (r])e~ z / 2 M 1 ^(2fce 2 ) ^ 
J B(w)c 1 (?7)e^/ 2 M_^ ^(2A;e z ) 



V 



-iB(v)c 2 (v)e~ z/2 M l ^r x (2ke z 



exp(ik x x), 



J 



where A{v) and -B(f) satisfy the system coupled system of equations 



d 



— -{k x -A 1 (y)))B = ikA, 



(81) 



(82) 



(83) 



± + (k x -A 1 (y)) S jA = ikB, (84) 

where v was defined in Eq. (|29|). 

The corresponding solution of Eq. (|59D in terms of the Whittaker functions Wk^(z) has 
the form 

^ / I^(w)c 1 (r/)e- 2 / 2 W r _ 1 /~(2£:e 2 ) ' { 



-iVA J B(u)ci(^)e- z / 2 W x /^(2A;e 2 ) 

~2> V A 

A(t»)c 2 (?7)e- 2 / 2 W r 1 /r(2A;e 2 ) 



V 



J B(w)c 2 (r/)e- 2 / 2 W r 1 



fee 2 



exp^/c^a;) 



(85) 



/ 



Let us look for solutions of the system Q8"3"D and (Q) when the electromagnetic potential has 
the simple functional dependence A±(y) = A±y. In this case one can obtain exact solutions 
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for A(v) and B(v) in terms of hypergeometric functions. After making the change of variable 



and using the recurrence relations [E? 



(b - l)M(a, b-l,z) = (b- l)M(a, b, z) + z 



dM(a,b,z) 
dz 



(86) 



1 dM(a, b, z) 
a dz 



+ M(a, b, z) = M(a + 1, b, z), 



(87) 



dU (a, b, z) 
dz 



— U (a, b, z) = —U (a, 6 + 1,2), 



we find that the general solution of the system of equations ( |83| ) and fl34| ) reads 

1 1 



.4 



/A- 



1 -1^2 

-e 2 



(88) 



(89) 



The exact solution of the system of equations 
interaction has the form 



CM-^r + llv 2 )). 



(90) 



AA X 2' 2' 
in the absence of electromagnetic 



(91) 



(92) 



fl = y/A: 2 fcg tfca- c ^f#Zj£ v _ yjk 2 k 2 x + %k x c ^ c -iy/&=ia v 
k k 



where C\ and C2 are arbitrary constants. The solutions of the Dirac equation ( |8"2"D and 
exhibit an asymptotic behavior which can be identified with the quasiclassical solutions of 
the Hamilton- Jacobi equation (§). With the help of the asymptotic expressions (|38|) , we 
find that the Dirac spinor $ as z — > —00, and y — ► 00. takes the form 



-UCi(77) 



e" ke "e^e"Tt)4 1 exp(zAvc) 



(93) 
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where the functions ci(ry) and 02(77) satisfy the system of equations ( j30|) and 
asymptotically large values of z we have that the spinor <3> takes the form. 



For 



zvAvci()7)( 



V 



e fee e 2 y2A! exp(z/c a ,a;) 



(94) 



-z2v^4Tc 2 (r7) 
-2kvc 2 (r]) J 

Looking at the solution of the Hamilton- Jacobi equation we can identify fl82|) and (j85|) as the 
corresponding quasiclassical modes as 2 — > — 00 and z — > 00, respectively. An approximate 
expression for the time dependence of the spinor $ can be obtained with the help of the 
WKB approximation. In this case we obtain 



c 2 (f?) ~ ci exp(iu;(77)), 



(95) 



ClO 



u(ri) + Ma (77) 



exp(ia;(77)), 



(96) 



where C\q is a normalization constant and uj{rf) = iM-^ + M 2 a 2 — A. Looking at (p5|)-(p6|) 



and (p3f) we readily see that, for large values of rj we obtain Ci{rf) — > — i 



ClQ 



2Ma(ij) 



exp(zu;(r7)). 



Analytic solutions of the system of equations ( p5|) and (|96| ) can be obtained for some par- 
ticular expansion parameter airf) P^JSH 



V. CONCLUDING REMARKS 

In this article, we have solved the Klein-Gordon and Dirac equations in an open cosmo- 
logical universe with partially horn topology. The solutions of the relativistic wave equations 
are expressed in terms of special functions. In Sec. IV we have shown that the algebraic 
method of separation p2|- p^j35[] permits one a complete separation of variables of the Dirac 
equation in the line element associated with a horn topology. The identification of the qua- 
siclassical modes with the help of the relativistic Hamilton- Jacobi equation shows that this 
method is a very useful tool in the study of quantum effects in curved spaces. 
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As a final remark, we should mention that the introduction of nonstandard topologies 
in order to describe the large scale structure of the space-time also opens new possibilities 
to discuss quantum effects in globally inhomogeneous and anisotropic backgrounds in the 
presence of non-trivial electromagnetic interactions. 
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